Transfer matrix eigenvectors 
of the Baxter-Bazhanov-Stroganov T2-model for N = 2 



Oleg Lisovyy 



in 
o 
o 

(N 
O 

Q 



School of Theoretical Physics, 
Dublin Institute for Advanced Studies, 
10 Burlington Road, Dublin 4, Ireland 

Bogolyubov Institute for Theoretical Physics 
14b Metrologichna str., 03143, Kyiv, Ukraine 



d 



> 

O 

in 
o 



> 

x 



olisovyyOstp . dias . ie 
Abstract 

We find a representation of the row-to-row transfer matrix of the Baxter-Bazhanov-Stroganov 
T2-model for N = 2 in terms of an integral over two commuting sets of grassmann variables. 
Using this representation, we explicitly calculate transfer matrix eigenvectors and normalize 
them. It is also shown how form factors of the model can be expressed in terms of determinants 
and inverses of certain Toeplitz matrices. 

1 Preliminaries 

The T 2 -model was originally introduced by Baxter in the work pQ, where it appeared in relation to 
the superintegrable case of the chiral Potts model. Later it was used by Bazhanov and Stroganov 
to establish a connection between six- vertex model and chiral Potts model [Hj. This connection has 
allowed to obtain a system of functional relations for transfer matrices of these models |2. and has 
led to the derivation of exact formulas for the free energy [3] and order parameter [I] of the chiral 
Potts model. Following the authors of [§1 Q2]> we will use instead of the name 'r 2 -model' the name 
'Baxter-Bazhanov-Stroganov model' (or simply 'BBS model'). 

BBS model is a system of spins, living on a square lattice and taking on N values 0, 1, . . . , N — 
1. The interactions exist only between nearest neighbours. In addition, the difference b 2 — 61 of 
neighbouring spins, living on the same vertical line (62 is higher than 61) is allowed to take on only 
the values and 1 (mod N). Consider an elementary plaquette of the lattice, drawn in the Fig. la. 
Boltzmann weights W(bi, b 2 , 63, 64), associated to this plaquette, are defined in the following table 
(our &!, 6 2 , b 3 , 64 correspond to d, a, b, c of and to 64, b x , b 2 , b 3 of 0): 
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a) b) 
Fig. 1. a) numeration of spins of an elementary plaquette b) graphical representation of the transfer matrix 



Here to = e 2 " K% l N and t,x,x',y,y', fi, fx' are parameters. It is easily seen that the model is Z^- 
symmetric: if one shifts the spins bi, . . . , 64 by 1, all plaquette Boltzmann weights remain unchanged. 

Consider now the case N = 2. Let us use instead of bi, ... ,64 new spin variables <7j = (— l) hj 
(j = 1, ... ,4), taking on the values ±1. The most general ^-symmetric Boltzmann weight is given 
by the following formula 



W(cri, cr 2 , cr 3 , 04) = a + ^2 a ' l i ° i(T 3 + fl4 Cr l°"20"30"4j • 



l<i<j<4 

The coefficients do, {a-ij}, Q>4, which correspond to BBS2 model, can be written as 



a = 


(y + fit)(y' + fi')/(4yy'), 




(y-fit)(y'-fi')/(Ayy'), 


QoQi2 = 


(y-fit)(y' + fi')/(4yy'), 


«0«34 = 


(y + lit)(y'-^)/(Ayy'), 


«0 a 13 = 


{l + xfi){t + x'fi')/{Ayy'), 


«0fl24 = 


(1 - x/i)(t-iV)/(%')- 


dQd\i = 


(l + xfj)(t-^fj!)/(4yy'), 


a a 23 = 


(l-x^)(t + x'fi r )/(4yy'). 



It was pointed out in [0] that these coefficients satisfy a 'free-fermion condition' 

O4 = Ol2 a 34 — a 13 a 24 + a 14 a 23- 



(1.2) 
(1.3) 
(1.4) 
(1.5) 
(1.6) 
(1.7) 
(1.8) 
(1.9) 
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Partition function of the model with plaquette weight (jl.lj) . satisfying the condition (JTTTDJ), was 
calculated by Bugrij [S] even in the case of a finite lattice. This result has allowed to obtain the 
eigenvalues of the BBS2 transfer matrix without solving any functional relations |9j. From the tech- 
nical point of view, the condition (jl.lOj) means that the Boltzmann weight (jl.lj) can be represented 
as the following integral over four auxiliary grassmann variables ip 1 , ip 2 , ip 3 , ip 4 : 



W(ai,a 2 ,(T 3 ,a i ) 



dip 1 dip 2 dip 3 dip 4 exp< auCi^ip 1 "^ 2 — a n oia- i ip l ip' i + ai i oio 4 ip 1 ip 4 + 



+a230"2cr 3 ^ 2 ^ 3 - a 2 4a 2 a 4 ip 2 ip 4 + a 3A a 3 a 4 ip 3 ip 4 \ e^WV*. 



(1.11) 



Throughout this paper, we will use the convention that 'dotted' grassmann variables commute with 
the usual ones, and that the variables inside each set anticommute: 

= ip^f 3 = -^> a , ip a ipP = ^ip a V a, (3. 

The method of grassmann variables was initially designed as a method of simple calculation of the 
partition function of the 2D Ising model. It was discovered and improved by different authors 
is probably the earliest reference). The use of two commuting sets of grassmann variables, which is 
crucial for our further discussion, was suggested in jjj. 

The main drawback of the method of grassmann integration is that it does not give the eigen- 
vectors of the transfer matrix, which are necessary ingredients in the computation of correlation 
functions and form factors. Even in the case of the Ising model, the only known practical way of 
calculation of these eigenvectors is the algebraic method of Kaufman ^1] (it should be mentioned, 
however, that recently a considerable progress has been achieved j!3| in the calculation of the eigen- 
vectors of the BBS transfer matrix using Sklyanin's method of separation of variables). It consists of 
two steps. First one should remark that the transfer matrix induces a rotation in a certain Clifford 
algebra. Then the eigenvectors are given by certain vectors from a Fock space, associated to the 
basis of this algebra, in which the above rotation is diagonal. Although Kaufman's method was later 
extended to some other free-fermion models [15 , it does not seem to work neither for the general 
free-fermion model nor in the case of the BBS2 model 1 . The main complication, as compared to the 
Ising model case, is that one should guess the explicit form of the appropriate rotation of the Clifford 
algebra. 

Having spent some time trying to guess the answer for the rotation, the author has finally found 
another method, which links grassmann integral approach with the transfer matrix formalism and 
allows to obtain the eigenvectors of the tranfer matrix of the general free-fermion model (i. e. the 
model with plaquette weight ()l.ip . satisfying the condition (|1.10Jl ). The present paper is devoted to 
the exposition of this method. 

This paper is organized as follows. In the next section, we find a convenient representation of the 
row-to-row transfer matrix of the periodic BBS2 model (or, rather, general free-fermion model) in 
terms of a grassmann integral, involving two commuting sets of variables (formulas (J2.13|) . ()2.15j) ). 
In Section 3, the eigenvectors of this transfer matrix are calculated (basic ansatz is given by (|3.3jl ). 
It should be pointed out that the form of the answer depends on whether the number of sites in one 
row of the lattice is even or odd. In Section 4, we find a dual basis of eigenvectors and normalize 
them. It is also shown that one can express form factors of the model in terms of determinants and 
inverses of certain Toeplitz matrices. Finally, in the last section the above results are specialized to 
two particular cases (BBS2 model and Ising model) and are rewritten in more common notation. 

1 Kaufman's method has also some drawbacks. First, it does not give a convenient representation of eigenvectors in 
terms of initial spin variables. Therefore, one is forced to do all the calculations in purely algebraic setting. Second, 
in the Ising case the transfer matrix spectrum is highly degenerate, so the eigenvectors are not determined uniquely. 
However, in the calculation of correlation functions and form factors one is typically interested in a very precise basis 
of eigenstates; in addition to the transfer matrix, they should also diagonalize the operator of discrete translations. 
Kaufman's method does not guarantee this last condition. 
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2 Grassmann integral representation 
for the transfer matrix 

Let us introduce the row-to-row transfer matrix of the BBS2 model. It is given by the product of 
plaquette Boltzmann weights over one row (see Fig. lb), 

L 

T[a,a'} = l[W(a J ,cr' j ,a' j+1 ,a 3+1 ), (2.1) 
i=i 

where periodic boundary conditions are imposed on spin variables: 

0"L+1 = 0"1, = 01- 

This matrix naturally acts in the 2 L -dimensional vector space V, composed of functions of L spin 
variables a%, . . . , vl- Namely, for any f[a] £ V we define the left action 

(Tf)[a] = Y,T[a,ayy]. 

Partition function of the BBS2 model onLxM lattice, wrapped on the torus, may be expressed in 
terms of the eigenvalues of T: 

Z(L, M) = Tr T M = £\ . . £ T[a^\ a®] . . . T[a^ M \ a^\. 

[<7«] [a( M )] 

In order to compute various correlation functions, one also needs to know matrix elements of local 
field operators in the normalized basis of eigenstates of T. To obtain these eigenstates in an explicit 
form, let us first find a convenient representation of the transfer matrix. 

From the formulas (jl.llj) and (j2.1j) it follows that one can write T in the form of a grassmann 
integral, 

T[cr, a'} = a o J V^Vip ex P ( fll2 a i a> i ~ fll3 a i a 'i+ 1 ^j^j+i + a ^ a j a j+i 

+a 23 ay j+1 ip]ipj +1 - a 24 Oj<r j+ i + ^34 Oj+i<rj+i tf+itf+i ) \ II e ^ e ^ e ^ +1 e ^ +1 ' 

/ J 3=1 

where the measure is given by 

L 

Vi{jViP = Y[ (d^) d^) dtf +1 
i=i 

Now let us make the change of integration variables: 

i)) -> o-j Vj, -> o-j 1%, ipj -> V>i °j V»i > j = 1, . . . , -L. 
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This change does not affect the measure, since every Oj and a'j appears in it twice. Spin variables 
disappear from the first (quadratic) exponential under the integral, but emerge in the 'tail'. Namely, 
one obtains 



T[a,a'\ = a^ VipVtp e 



,Si[V,V>] 



i=2 



(2.2) 



where Si[ip,ip] can be schematically represented as 



;2.3) 



with ( ip up J = ( ip l ip 2 ip 3 ?/> 4 ) and 



( a 12 -a 13 V x a 14 V x \ 

o 12 a 23 V x -a 24 V x 

-au^-x -a23^7-x a 34 

\ -a 1A V- x -a 24 V_ a; a 34 / 



(2.4) 



Here V x denotes the operator, shifting the lower indices of grassmann variables by 1. It obeys 
periodic boundary condition (V X ) L = 1. For example, one has 

L-l 

Note that in (|2.2j) we have rearranged the tail, assembling together the exponentials, containing 
the same spin variables. It seems, however, that the exponentials e ' 1 ^ and e CT1 ^i are not on their 
'right' places. One may correct this, observing that for any function F[ip] and for any grassmann 
variable ip a we have the identity 



F[if)] e* a = 
Now, introducing the notation 



2 2 



F 1 [il>] = e°^ J] 
i=2 



i=2 



one may rewrite the tail as 



x< e 



F 2 [-j>] 



;2.5) 



Expanding this last expression, one obtains 16 terms. However, some of these terms are equivalent, 
since the simultaneous change of the signs of all ip and ip does not affect the value of the integral. 
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Then one may easily check that ()2.5|) may be replaced (after appropriate change of variables) by the 
following combination, containing only 4 terms: 

F x \fj}\ e CTl ^F 2 [^] -> X - {e^F^] e^F^] + e^FiHA] F 2 [V>]+ (2.6) 

+e CTl ^F 1 [V>] e^i^] - e- CTl ^Fi[-V>] e^F^]} . 

The third term of the last expression has the desired form and there is no need to transform it 
further. If we make the substitution ip\ — > —ipf, ipf — » —ipf in the integral, corresponding to the first 
term, it will have almost the same structure. The only difference is that the boundary condition for 
the shift operator V x becomes antiperiodic: (V X ) L = —1. 

Next one should remark that the second and the fourth term in (|2.6|) can be obtained from the 
first and the third one, respectively, by changing the signs of the spins 01, . . . , ol- This change can 
be realized, using the operator of spin reflection U, whose defining property is that {U f)[a] = /[— a] 
for any vector f[a] G V. Matrix elements of U may be explicitly written as 

3=1 

Summarizing the above observations, we obtain the following representation for the transfer 
matrix: 

T=I±^T JV5 +— -T R , (2.7) 
2 2' K J 

where 

T NS(R) a>] =a L f v ^ v ^ exp | s ns(r) ^ ^ | ^ ^>,-> ( »,^ Y[ (e^e°^) , (2.8) 

J 3=1 3=1 

and both actions S^ 8 ^ [ip, ip] are defined by the formulas ()2.3|) - (j2.4|) . Upper indices NS and R 
correspond to antiperiodic (Neveu-Schwartz) and periodic (Ramond) boundary conditions, satisfied 
by the shift operator V x . 

The matrices P± = have the properties of projectors, i. e. P± = P±; thus their eigenvalues 
are equal to either or 1. The eigenvectors, corresponding to zero eigenvalues of P + (P-), are odd 
(even) under spin reflection. It means that a vector f[a] G V will satisfy (P + f)[a] = ((P_/)[cr] = 0) 
iff f[a) = —f[—a] (respectively, f[a] = f[—a]). Analogously, the eigenvectors of P + (P_) with 
eigenvalue 1 are even (odd) under spin reflection. 

The operator U commutes with the transfer matrix T. Therefore, these two matrices can be 
diagonalized simultaneously and one may choose the eigenvectors of T so that they are either even 
or odd under the action of U. Let us take an even eigenvector f e of T, and denote by A/ e the 
corresponding eigenvalue. Acting on f e by both sides of the relation (|2.7jl . and using the fact that U 
commutes with T NS and T R as well, one obtains 

Xfje = Tf e = (T NS P + + T R P_)f e = T NS f e , 

that is, any even eigenvector of T is an eigenvector of T NS with the same eigenvalue. Similarly, any 
odd eigenvector of T is an eigenvector of T R . Conversely, any even eigenvector of T NS and any odd 
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eigenvector of T R are eigenvectors of T. Therefore, the set of all transfer matrix eigenstates splits 
into two parts: NS-sector (even eigenvectors of T NS ) and R-sector (odd eigenvectors of T R ). The 
problem of diagonalization of T is then reduced to the calculation of eigenvectors and eigenvalues of 
matrices T NS and T R , given by the formula (|2.8|) . 

Remark. Having diagonalized T NS and T R , one also gets for free the solution of the BBS2 model 
with antiperiodic boundary conditions for spin variables (in one direction). It is easy to understand 
that the transfer matrix of such model is given by 

rpa _ 1 ~ ^ rpNS , ^ + U rpR 

2 2 

and, therefore, the set of its eigenstates is composed of odd eigenvectors of T NS and even eigenvectors 
oiT R . 

The representation (J2.8)) can be simplified even further by integrating over fermionic degrees of 
freedom, which are not coupled to spin variables. We mean the following: elementary factors from 
the products of ()2.8|) can be written as 

Let us now introduce instead of if) and ip new grassmann variables 

Since the jacobian of the transformation ()2.9|) is equal to 1, the integration measure transforms as 

L 

Dip Dtp — > Dip T>(p T>7] T>r] = J^J (d(fj d(pj drjj drjj j . 
Then the integral (J2.8|) may be rewritten as 



rpNS(R) j- 



a, a = a 



r L L 

/ vcpVifVvvr] exp{^ 5w [^ <p, 77,77]} ( 2 - 10 ) 



where the action S^ 5 ^^, 0, 77, 77] is given by 

S? S{R) [<p,<p,T],ri] = 2 ( ^ ^ ^ v)D 2 (f <p V v) T , 

( ai2 -ai2-ai3Va; -l + auVx \ 

012 -l + a23V x -ai2-a24V x 

-ai2-ai3V_ x l-fl23V- x -a 2 3(V :c - V_ x ) a i2 + a 34 + ai 3 V _ x + a 2 4V 2 

\ 1 - au V_ x -ai 2 - a 2 4 V_ x a i2 + a 34 + a^V^ + a 2 ^ -x -au(V x - V_ x ) / 

Let us now integrate over 77 and rj in the representation (|2.10jl . This integration can be done relatively 
easily, since , (p, 77, 77] is diagonalized by Fourier transformation. Namely, if one denotes 

1 L 

( ¥p ¥p V P V P ) = -?f^2 e ~ WJ ( fi Vj Vj ) , ( 2 - n ) 

V j = i 
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D 
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X 



then 

S2 S{R) [<p,<j>,V,v] = g Yl ( } ( V-p <P- P V-p V-p ) D 2 (p) ((p p P % V P ) T , (2-12) 

p 

where the one-mode matrix D^p) is given by 

/ a 12 -ai2-«-i3e ip -l + ai 4 e ip \ 

f) ( \ _ °12 -1 + a 2 3 e ip — ai2 - a 24 e ip 

— ai2 - a 13 e" ip 1 - a 23 e~ ip —2i a 23 ship a i2 + a 34 + ai 3 e~ ip + a 24 e ip 

V 1 — ai4 e~* p — a 12 — a 24 e~ ip a 12 + a 34 + ai 3 e ip + a 24 e~ ip — 2i a 44 ship y 

and the indices NS and R in the sum (j2.12|) mean that the corresponding quasimomenta run over 
Neveu- Schwartz values (p = (j + ¥)i 3 = 0, 1, ... ,L — 1) or, correspondingly, Ramond values 
(j) — j, j = 0, 1, . . . , L — 1). Note also that the integration measure can be written as 

nNS(R) / \ 
( difp difp drip drjp J 

p 

Thus the 2L-fold integral over rj and rj in the representation ()2.1()j) factorizes into a product of 
4-fold (over r)± p , fj± p ) and 2-fold integrals. Double integrals correspond to the mode p = (always 
present in the Ramond sector) and p = n (present in the Ramond sector for even L and in the 
Neveu- Schwartz sector for odd Li). After a little bit cumbersome but nevertheless straightforward 
calculation one then obtains 

T NS W[a,a'] = ( NS ^ I V NS ^V N3 ^ exp{S NS ^,^]} f[e^ f[e^ } (2.13) 

** „■ i „■ 1 



where 



i=i j=i 

v ns ^ 9 = i[ NS(R) dipp , = i[ NS{R) dv P , 

p p 

C NS ^ = a L \[ NS{R \y\ (2.14) 

p 

X P = W12 + «34 + (oia + «24) cosp] 2 + [(a 13 - a 24 ) 2 + 4ai 4 a 23 ] sin 2 p, 
and the action S NS ^ [cp, 0] is given by 

(2 ' 15) 

with 

XpGn(p) = 2isinp a 23 + a 12 a 24 + ai 3 a 34 + a 14 a 4 - 2(a 14 a 23 - a 13 a 24 ) cosp 
Xp G22(p) = 2i ship a i4 + ai 2 ai3 + 0-34024 + a-2304 — 2(ai 4 a 2 3 — 013O24) cosp 
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X P G 12 (p) = X P G 2 i{-p) 



+ 
+ 



(oi 2 + a 34 )(a 4 + 1) - (a i4 + a 23 )(ai 3 + o 24 ) 
(ai 3 + 024) (04 + 1) - (ai4 + 023) (012 + o 34 ) 
(a 24 - ai 3 )(a 4 - 1) + (o 44 - o 23 )(oi 2 - 034) 



+ 

cosp + 
i ship . 



As we will see in the next section, this final representation for T NS ( R > (given by the formulas (|2.13jl . 
(|2.15]1 ) allows to obtain all transfer matrix eigenvectors almost immediately. Concrete form of the 
functions = 1,2) does not play any essential role. 

For further convenience and making parallels with the work [Oj, let us also introduce the notation 
v P = 4:X P G 12 (p) and 



COSJ9 . 



up = 2 X J 1 - G n (p)G 22 (p) + G 12 (p)G 21 (p) 



;i+a 4 ) 2 +(a 12 +a 3 4) 2 +(ai3+«24) 2 +(ai4+a23) 



+4 



(^12+034) (013+024) - (014+023) (I+04) 



Remark. It should be pointed out that we did not care about the correct overall sign of T NS ^ in 
the representation (|2.13|) . However, using the fact that for a i2 = a 13 = au = o 2 3 = a 24 = a 34 = all 
the eigenvalues of T should be equal to 2 l Qq , one can restore this sign at any stage. 



3 Transfer matrix eigenvectors 

It appears that the form of the eigenvectors of T NS and T R depends on whether L is even or odd. 
Moreover, quasiparticle interpretation of the eigenvectors and eigenvalues is different in different 
regions of parameters of the BBS2 model. Below we will consider various cases in order of increasing 
difficulty. 

3.1 NS-sector, even L 

If L is even, then the Neveu-Schwartz spectrum of quasimomenta does not contain the values and 
7r (the only values with the property p = —p mod 2tt). 

The simplest ansatz for an eigenvector f[a) G V of the matrix T NS is given by an integral over L 
auxiliary grassmann variables £1, . . . , £1,: 

V NS i exp |£ 2 £_ p A(p) £ p | j[ < (3.1) 

Here A(p) = — A(— p) is an unknown odd function to be determined, and {£ p } denote Fourier 
components of £. The indices 4p and ^ in sums and products will be used to indicate that the 
corresponding operations involve only those Neveu-Schwartz and Ramond quasimomenta, which lie 
in the open interval (0, n) (for the NS-sector and even L, this is exactly one half of the Brillouin 
zone). Note that for even L the function (|3.1j) is even with respect to the action of U: the reversal 
of all spins is equivalent to the change of variables £ — > — £. 
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Let us now act on f[a] by the matrix T NS . Since the fields <p and £ commute, the sum 
T NS [cr, a'} /[a'] can be easily evaluated and one obtains 



(T NS f)[a] = 2 L ( NS j V NS ipV NS <p e sNS[ ^ ] f[e^ x 

3=1 



X 



/( NS_ L \ 

V NS i expj^ ^ A (?) & + E & ^ ] • 



p j= 
After integration over £ one finds the exponential of a quadratic form in (p, 

L 



/ NS ^ 
V NS i exp IJ2 2 & + E & & 

I p i=i 

j V NS £ exp ^(p) £ p + 0_ p £ p + P £- ? 



NS \ NS 



] [ 2 4(p) exp <^ 2 (p_ v A 1 (p)if p 



which can then be pulled through the 'linear' exponentials. Then one may integrate over (p and 
obtain 



L 



ns , . r \ EM. 

(T NS f)[a}=2 L C s H 2 (l-A(p)G 22 (p)) /P^exp £ 2 ^ p A'(p)^ Ijje'. 

p I p J i=i 

with 

-4 (p) = Gn 0) + — , . ■ 

1 - A(jo)G 22 (p) 

Therefore, the function (j3.1|) will be an eigenvector of T NS iff for all NS- values of p from the interval 
(0,7r) one has A(p) = A'(p). This equation is quadratic in A(p), and its roots are given by 



^ 1 + 6^)^(1-0^-40^) 
A (P) " 2~GM ' (3 - 2) 

where we have introduced the notation 

G(p) = G 11 (p)G 22 (p) - G 12 {p)G 21 (p) . 

Thus the formula ()3.1|) gives 2 i//2 eigenvectors of T NS , corresponding to different choices of the set 
of one-mode roots. 

One can take, for instance, A(p) = A + (p) for all p G (0,7r). The vector, corresponding to 
this particular choice, will be denoted by \vo,c)nsi since under some conditions, satisfied by the 
parameters of the BBS 2 model, it corresponds to the eigenvalue with maximum modulus. Similarly, 
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if we choose A(p) = A~{p) for some values p±, . . . ,pk G (0, ir), and A{p) = A + (p) for all the other 
NS-quasimomenta from the interval (0,7r), then the corresponding eigenvector will be denoted by 
\pi, —pi] ■ ■ - Pk, —Pk)NS- The origin of this notation will become clear soon. 

In order to find all the eigenvectors of T NS , only a slight generalization of the ansatz (|3.1|) is 
needed. Namely, let us define 

jvs L 

V NS t W F^&lle'* (3.3) 

p 3=1 

Here each of the indices {i p } can take any of the four values, which we will conventionally denote by 
1, 2, 3, and 4. Corresponding functions -Fj(£_ p ,£ p ) are defined as follows: 

*i = exp^ p A + (p)^, (3.4) 

F 2 (Z- P ,Q = £- P , (3.5) 
Fz(^ p ,Q = f p , (3-6) 

*4(£-p,&) = exp^_pA-(p)^). (3.7) 

Similarly to the above, one should act on /{^M by T NS , then to sum over the intermediate spin 
variables, to integrate the result over £ and, finally, over ip. Then it is straightforward to verify that 
the formulas (J3.3|) - ()3.7j) indeed define an eigenvector of T NS with the eigenvalue 



JVS 



A NS } = 2 L C NSY[> (3 . 



P 



where 'one-mode' eigenvalues are given by 

X 1 (p) = l-A + (p)G 22 (p), (3.9) 
A 2 (p) = G 12 (p), (3.10) 

Hp) = G 21 ( P ), (3.ii) 

A 4 (p) = l-A-(p)G 22 (p). (3.12) 
One can also rewrite them in the following way (see the end of the previous section for the notations): 



Ai(p) = ^- , x *(p) = ^ x ^ = ^ A4 ^ = • (3 ' 13) 

The total number of found eigenstates is equal to 4 L//2 = 2 L , and thus the diagonalization of the 
matrix T NS is completed. 

Let us now turn to quasiparticle interpretation of eigenvalues and eigenvectors. It follows from 
(|2.14j) . (|3.8)L ()3.13|) that the eigenvalues can be written in the form 



where 

r-NS 



A ^ax = a oYl PP 2 ' Pp = U p + Ju 2 p -V p V_ p . (3.15) 



P 



Remark. The expression — v p v_ p is a quadratic polynomial in cos p. For the sake of simplicity, 
it will be assumed that the parameters { a ij} 1 <j < j< 4 are all real and chosen so that this polynomial 
has no roots inside the interval (—1, 1) (for example, this condition is satisfied, if one takes a%2 = 034 
and ai3 = 024). It means, in particular, that u p — v p v_ p is non-negative and has local extrema only 
at the points p = and p = 71. 

Consider also the operator of translations in discrete space R. Its action on an arbitrary vector 
f[a] G V is defined as 

(Rf)(a 1 ,a 2} ...,a L ) = f(a 2 ,a 3 ,...,a l ). (3.16) 

This operator commutes with the transfer matrix T, with the matrices T NS and T R , and also with the 
operator U of spin reflection. Since we have already diagonalized T NS and obtained nondegenerate 
spectrum, the eigenvectors ()3.3j) should diagonalize R as well. Actually, it is not difficult to verify 
that 

NS NS 

wg})M= n 2 eP n 2 e_ip /(KM- 

p i p =2 p\i p =3 

Now it is clear that the eigenvectors of T NS can be labelled by the collections of distinct NS- 
quasimomenta and interpreted as multiparticle states. One-particle energy is given by 

e(p) = — In — . 

Pp 

It may have a non-zero imaginary part, which is a general consequence of the fact that the transfer 
matrix T of the BBS2 model is not symmetric. The eigenstate, which contains particles with the 
momenta pi,...,p k , will be denoted by \pi, ■ ■ ■ ,Pk)NS- in order to determine, which one of the 
functions (|3.3|) gives the explicit form of this vector, one should decompose the set of momenta of 
particles from the state \pi, . . . ,Pk)NS into three parts: pairs of the form ±Pj, 'unpaired' momenta 
from the interval (0, n), and 'unpaired' momenta from the interval (it, 2tt). Then in the ansatz ()3.3|) 
one should set 

• i p — 4, if ±p appears in the first part, 

• i p = 2, if p appears in the second part, 

• i p = 3, if — p appears in the third part, 

• i p = 1 for all the other values of p. 

This procedure establishes the correspondence between the formulas (|3.3|) and usual quasiparticle 
notation. 

Remark. Recall that only even eigenvectors of T NS diagonalize the full transfer matrix T as well. 
It means that the total number of appearances of i p = 2 and i p = 3 in the functions ()3.3|) should be 
even. In other words, NS-eigenstates of T should contain even number of particles. 
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3.2 NS-sector, odd L 

When L is odd, the Neveu-Schwartz spectrum of quasimomenta contains the value p = it. To take 
into account this special mode, it is sufficient to slightly modify the ansatz (|3.3|) . Let us consider 

/ NS L 
v Ns i n 2 f ^-p,q n ^ ( 3 - 17 ) 

V 3=1 

where all the indices {i p }, except i w , take on four values as above, and the functions F\ . . . F4 are 
given by (|3.4)l - (j3.7jl . The index i n can have only two values, 1 and 2, and the corresponding functions 
Fx and F2 are simply 

Fx(o = i, A(o = e- 

One may verify that the function ()3.17j) gives an eigenvector of T NS with the eigenvalue 

A^ } = 2 i C JV5 ^( 7 r)n 2 A^Cp), 

p 

where all Aj(p) are defined as above and 

Ai(tt) = 1, A 2 (7r) = G 12 (7r) = -^. 

Since the total number of eigenvectors (|3.17jl is equal to 2 x A^ L ~ 1 ^ 2 = 2 L , the diagonalization of 
T is completed. 

The first thing that may seem unusual is that if we set i p = 1 for all p, including p = 71, 
the corresponding eigenvector of T NS will not always represent the physical vacuum. Moreover, this 
vector is odd under spin reflection and, therefore, it is not an eigenvector of the full transfer matrix T. 
Note also that 

P y 2 = 2 x y 2 max{M ij. 

Therefore, if one tries to write the eigenvalues in the form, analogous to (|3.14jl . then the result will 
be different in different regions of parameters. Namely, for 1^1/4%^ > 1 one obtains 



2 — i n jvS NS NS 



AS=As e ^ n 2 n 2 n n 



2 ^p ^— p 
2~ ' 



pfe=2 P P p\i p =3 Pp p\i p =4 P P 



and for \v k \/Axti < 1 we have 



A TVS _ A TVS / U T \ TT 2 TT 2 ^-p TT 2 Up U-p 

H P } ~ A — I — J 11 — 11 — 11 —2—, 

p|i p =2 ^ p|i p =3 ^ p|i p =4 P P 

where A^f x is defined by the formula (|3.15|) . Thus one can again interpret the eigenvectors of T NS 
as multiparticle states \p\, . . . ,Pk)NS- The main differences with the previous case are the following: 

• If |f 7r |/4x 7r > 1, then the states, containing a particle with the momentum p — tt, are given by 
the ansatz (|3.17|) with i n = 1; for 1^1/4%^ < 1 they correspond to the choice v = 2. 

• For 1^1/4%^ > 1 the eigenstates, which are even (odd) under spin reflection, contain even (odd) 
number of particles, while for |t> 7r |/4x 7r < 1 this number should be odd (even). 
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3.3 R-sector, odd L 

The treatment of this case is completely analogous to the previous one, since for odd L Ramond 
spectrum contains only one 'special' mode p = 0. All the eigenvectors and eigenvalues of the matrix 
T R are given by 

/&>m = / vR z ^0(6) n f w-p'tp) n e< ^> ( 3 - 18 ) 
j P j=i 

p 

Here the indices {i p } p ^o take on four values, i = 1,2, the functions {Fj}, {Fj}, {Xj} are defined as 
above and 

A!(0) = 1, A 2 (0) = G 12 (0) = -^_. 

Again, since we have 

1/2 1/2 J I vol 
Po = 2 Xo max ' 



Axo' 

the quasiparticle interpretation of eigenvalues and eigenvectors is different in the regions |vo|/4xo > 1 
and |vo|/4xo < 1- Namely, one has 

^ - (^T n ? ? n f ^ n f ^ *» w/** > i, 

where the eigenvalue with the maximum modulus, A R ax , is given by 

A* =ttV /2 - 

max J_ J_ rp 

P 

Similarly to the above, let us denote by \p±, . . . ,Pk)R the eigenstate of T R , containing k particles with 
distinct R-momenta Pi, ■ ■ ■ ,Pk- 

Note that even (odd) eigenstates of T R should contain even (odd) number of particles for 
|vo|/4xo > 1; an d odd (even) number of particles for |vo|/4xo < 1- This change can be easily un- 
derstood if we take, say, |vo|/4xo < 1) an d consider two eigenstates, \pi, . . . ,Pk)n and \0,pi, ■ ■ ■ ,Pk)n 
(pj ^ 0, j — 1, . . . , k). Then let us gradually increase |f |/4xo- When this parameter approaches the 
critical value 1, two eigenstates correspond to the same eigenvalue, and when it exceeds 1, the roles of 
two vectors swap around: the particle with zero momentum disappears from the second vector (thus 
decreasing the number of particles by 1) and appears in the first (the number of particles increases 
byl). 
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3.4 R- sect or, even L 

Since for even L the Ramond spectrum of quasimomenta contains both p = and p = it, the 
eigenvectors and eigenvalues of T R in this case can be written in the following way: 



J P j=i 



(3.19) 



Ai } = 2 L c R \ l0 mMii 2 \ l M- 

v 

From the physical point of view, here one should distinguish four different regions in the space of 
parameters. They have the following properties: 

• | Vol > 4xo 5 |^tt| > ^Xtt- The eigenstates of T R , containing a particle with the momentum p = 



(p = n), are given by the formula ()3.19|) with iq = 1 (v = 1). The eigenvectors, which are even 
(odd) under spin reflection, should contain even (odd) number of particles. 

I^o | > 4xo, I^ttI < ^Xtt- The eigenstates, containing a particle with the momentum p = (p = 



7r), correspond to io = 1 (k = 2). Even (odd) eigenvectors contain odd (even) number of 
particles. 

I Vo| < 4xq, \v w \ > Particle with the momentum p = (p = n) corresponds to io = 2 



(i n = 1). Even (odd) eigenvectors contain odd (even) number of particles. 

• |i>q| < 4xoj I^ttI < 4xtt- Particle with the momentum p = (p = n) corresponds to i = 2 
(in = 2). Even (odd) eigenvectors contain even (odd) number of particles. 

4 Norms and form factors 

In the present section, the problem of computation of correlation functions of the BBS2 model is 
addressed. Local fields will be represented by spin variables (i = 1, . . . , L; j — 1, . . . , M). In the 
transfer matrix formalism, 2fc-point correlation functions (o , i 1 j 1 o , i 3i j 3 . . . &i 2k ,j 2k ) can be written in the 
following way 2 : 

(cr iujl ...(T i2k j 2k ) = (4.1) 
= Z-\L,M) a h ^' 2 ~ J VV (2) ] a { £ T^VV (3) ] . ..a™ T 1 

[crM] [a< 2fc )] 

where it was assumed that j± < ji < ■ ■ -3ik- Let us introduce spin operator 



^M-0 2fc -ji)^(2fe) )(T (l)j ^ 



acting on functions f[a] G V from the left in the usual way. If we make use of the translation operator 
R (see formula (JHHEJ)) to define 

S itj = Ti- 1 R 1 - 1 3 lt i R 1 ' 1 T 1 "- 7 , i = 1, . . . , L, j = 1, . . . , M, 



2 All (2k + l)-point correlation functions vanish due to Z2-symmetry of the model. 
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then one may rewrite (j4.1j) as 

TV ( q. . q. . q. . r pM\ 

I V _ 11 \°ll,n°l2,j2 ■ ■ ■ °l2k,j2k 1 ) (a r>\ 

\ a h,h ■ ■ ■ a i2k,hkl — ^rpM ' V / 

Since all the eigenvalues of the transfer matrix T are known, the problem reduces to the calculation 
of the trace in the numerator. One would want to compute this trace in the basis of eigenstates of 
T. However, such computation is not quite straightforward, since the transfer matrix of the BBS2 
model is not symmetric and thus its eigenvectors are not necessarily orthogonal. Therefore, in order 
to construct the dual basis, one should separately find the eigenvectors for the right action of T, 
(since T is not symmetric, they can not be obtained from the eigenvectors, found in the previous 
section, by simple transposition). 

Assume for a moment that we have found all 'left' and 'right' eigenvectors of T. Let us denote 
them by \n) and (n\ , where n is any convenient set of quantum numbers, identifying the eigenstate 
(for example, the number of particles and their quasimomenta) . The resolution of the identity matrix 
in this basis of eigenstates has the form 



n 



b n \n)(n\, b n = l/(n\n). (4.3) 



The relation (|4.3|) means, in particular, that the trace of any matrix X can be written as 

TrX = ^b n (n\X\n) . 

n 

Recall also that the eigenvectors of T diagonalize as well the translation operator R. Therefore, 
inserting the resolution of the identity matrix into the representation (|4.2|) k times, one can rewrite 
2/c-point correlation function in the form of the so-called form factor expansion. For example, for 
the 2-point correlation function one has 



b m b n {n\S lt i\m)(m\S ltl \n) e-^^-J'O-^^-ia+iO-H^m-PnJCia-ii) 



( cr i 1 J 1 C r i 2 ,j 2 ) — „—M E n ' (4-4) 



Matrix elements {n\S\j\m), entering this formula, hereinafter will be referred to as form factors. 
Parameters E n and P n have the meaning of energy and total momentum of the state \n) (and (n|). 
They are related to the eigenvalues of T and R in the following way: 

T\n) = A max e- En \n)i R\n) = e iPn \n), 

where A max denotes the eigenvalue of T with the maximum modulus. In the BBS2 model, the values 
of E n and P n , corresponding to the multiparticle state |n) = \p±, . . . ,Pk), are given by the sums of 
one-particle energies and momenta. 

The generalization of the form factor expansion (j4.4j) to the multipoint case is straightforward. 
Thus in order to find all correlation functions, only three further steps should be made. First one 
should find the eigenvectors for the right action of the transfer matrix T, i. e. the functions f[a] G V 
such that f[ <7 ]T[o', a'] = A//[cr']. Then one needs to compute scalar products (n\n). Finally, the 
H 

most difficult task is the calculation of form factors (nlSi^lm). All these problems are treated (the 
third one with only a partial success) in the following subsections. 
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4.1 Eigenvectors for the right action of T 

The variables [a] and [a'} enter into the representation (j2.13|) for the matrices T NS and T R in a 
similar way. Therefore, one may construct the eigenvectors for the right action of T along the lines 
of Section 3. However, there exists even more straightforward way to obtain them. Note that the 
right action of T on f[a] G V coincides with the left action of the transfer matrix T = T T of another 
BBS2 model (see Fig. lb), characterized by the parameters 

(2 12 = ai2, ai3 = 024, «14 = <223) ^23 = a 14; «24 = a 23) &34 = &34- (4-5) 

Thus the eigenvectors we are looking for may be obtained from already found ones by the substitution 
(|4.5j) and mathching the eigenvalues. 

It is easy to verify that under the above substitution various quantities, used in the construction 
of eigenvectors and eigenvalues, change as follows: 

Gn(p) G l2 (p) \ ( G u (p) G 12 (p) \ ( G 22 (p) G 21 (p) 
G 21 ( P ) G 22 (p) J \ G 21 (p) G 22 (p) J V G M G u (p) 



l + G(p)T^(l-G(p)Y-AG 12 (p)GM 
A( P )^A (p) = ^-^ , (4.6) 

Xp * Xp Xp > * ) * ^p ^— p • 

Let us now consider, for instance, NS-sector and assume that L is even. Let fjigyW] denote the 
'left' eigenvector (|3.3j) . corresponding to the eigenvalue Under the substitution (|4.5j) 'partial' 

eigenvalues A 2 (p) and A 3 (p) (formulas (|3.10j) . (|3.11jl ) exchange their roles, while \i(p) and A 4 (p) 

r {i P }l 

to the same eigenvalue as /^[o - ], is given by 



remain unchanged. Then it becomes clear that the 'right' eigenvector f$^[o] of T NS , corresponding 

eNS I 



/ jVS 
p i=i 

with 

^i(tp,4) = exp(t P i+(p)4), (4.8) 

^ 2 (e-p,4) = 4, ( 4 -9) 
F 3 (Lp,i P ) = i- P , (4.io) 

MLpAp) = exp^_ p i-(p)4), (4.11) 

the functions ^(p) being defined by the formula (|4.6|) . Dotted grassmann variables £ are used in the 
representation ()4.7|) for further convenience in the computation of scalar products and form factors. 

In order to obtain a similar answer for the other cases (Neveu-Schwartz sector for odd L and 
Ramond sector), it is sufficient to substitute in ()3.17|) . (|3.18|) and ()3.19|) instead of F1...F4 new 
functions Fx . . . F 4 . The functions Fi and F 2 , which are responsible for the special modes p = 0, tc, 
remain unchanged. 
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4.2 Normalization 



It is instructive to consider not only the norms (n\n), but also general scalar products (m\n), and to 
verify by hand that (m\n) = for m ^ n. First one should remark that the eigenvectors of T, which 
belong to different sectors, are orthogonal, since they correspond to different eigenvalues of U . Thus 
one may look at each sector separately. Let us now consider, for instance, the Neveu-Schwartz sector 
for even L. Let us take a 'right' eigenvector ^({zp}) d = f{ip}l a ] (gi ven by the formula (j4.7|) ) and 
a 'left' eigenvector \{j P }) NS = f^i* 7 } (given by the formula (jSISj) ), and then compute their scalar 



product 

NS { 



Since the fields £ and £ in the representations ()3.3|) and ()4.7j) commute, the summation over inter- 
mediate spins can be easily done and one obtains 



/NS , . v 

p 

Calculation of this factorized integral gives 

NS 

W s({ ? p}|{jp})jvs = 2 L II 2 a w(P)' 

p 

where the functions cxij{p) can be assembled into a 4 x 4 matrix 



/ A + (p)A + (p) - 1 A + (p)A-(p) - 1 \ 

-10 

0-1 

\ A~(p)A + (p) - 1 A-(p)A-{p) - 1 / 



Using the explicit formulas for A (p) and A ± (p), one may check that A ± (p)A Zf (p) = 1. Therefore, 
'right' and 'left' eigenvectors, corresponding to different eigenvalues, are orthogonal (as it should be). 
The norm NS ({i P }\{i P }) NS is given by 



NS 

N s({h}\{h}) NS = 2L U 2 <*iM ( 4 - 12 ) 



where we have introduced the notation 

ai(p) = 1 - A + (p)A + (p), a 2 (p) = a 3 (p) = 1, a 4 (p) = 1 - A"(p)A-(p), 

and corrected the overall sign. One may check that the answer for the Neveu-Schwartz sector and 
odd L is given by the same formula (|4.12jl . The only things that change in the Ramond sector are 
the values of quasimomenta. 
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4.3 Form factors 



Since the eigenstates of T from the same sector are all simultaneously even or odd under spin 
reflection, all form factors of type NS-NS and R R are equal to zero. Now assume for defmiteness 
that L is even and consider a 'right' eigenvector ^({^p}] — /{I'}! "] from the Neveu-Schwartz sector 
and a 'left' eigenvector [{jp})^ = f{j p }[ a ] fro m the Ramond sector. Let us calculate the form factor 



NS 



R 



After summation over intermediate spins one obtains 

NS ({h}\Si,i\{j P }) R 



ii) exp 



(4.13) 



. k=l 



/R NS 
p i 

Unfortunately, we have not managed to find a compact expression for this gaussian integral, 
although we strongly suspect it is possible. In order to illustrate emerging difficulties, let us assume 
that | z; | < 4xo, I IV I > ^Xn (this region mimics ferromagnetic phase), and consider the simplest 
possible form factor Ns{vac\Si t i\vac) r, which corresponds to the following choice: jo = 1, jn = 2, 
j p = 1 for all p G (0, 71"), i q — 1 for all q G (0, it). One then obtains 



Ns( vac \Si,i\vac) 



R 



/( R JVS L j 

I P q fc=l J 

Quadratic form in the exponential consists of three pieces, which can not be diagonalized simulta- 
neously: the first and the second piece are diagonal in the Fourier basis with Ramond and Neveu- 
Schwartz values of discrete quasimomenta, and the third one is diagonal in the coordinate repre- 
sentation. Actually, one can now remove the dots, using the following rule: all quadratic terms in 
the exponential, containing a dotted variable on the left, should change their signs. Performing this 
operation and passing to the coordinate representation in all terms, one obtains 



Ns(vac\Si^\vac) 



R 



^JvtVv J2(-l)% (6 + m) exp ji ( £ V ) 



A + 1 
-1 -i H 



where antisymmetric L x L matrices A + , A + are given by 

A t*> = \H Ra+ (p) ^~ x ' ] > A **> = jY, NSa+ ^ el9ix ~ x ' ] > 

p^0,7r q 

Evaluation of the gaussian integral (|4.14jl gives 

2 L 



(4.14) 



x, x 



L. 



Ns{vac\Si : i\vac) 



R 



Pf (i + ) Pf (H) J2 



k=l 



H- 1 



A" 



H- 1 



(4.15) 



ik 
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where L x L matrix H is also antisymmetric and has Toeplitz form: 

H xx , = - A + xx , = \ Y? S A-(q) e*'-* - \ £ V(p) 1 

Thus the problem of computation of the form factor NS (vac\Si j i\vac) ^ is reduced to the calculation 
of the determinant of H and inverse matrix if -1 . Actually, this is also the case for more complicated 
form factors. In spite of the remarkably simple form of the matrix H, we have not succeded in the 
calculation of Pf(H) and H~ x . However, we believe that the representations of type (|4.15|) are still 
useful, since they effectively reduce initial 2 L -dimensional problem to an L-dimensional one. 

It should also be pointed out that form factor Ns{vac\Si t i\vac) r enters into the definition of the 
order parameter of the BBS2 model. More precisely, one has 

z n2 def r ( v 1 A r N s(vac\S 1A \vac) RR (vac\S 1A \vac) NS 
(cr) = hm hm {cr hl a Lj } = hm ■ — '— . 

ij--Kx> \l,m-+<x> ' J L^tx> Ns{vac\vac)Ns R\vac\vac)R 

Although we have not managed to obtain a closed expression for this form factor, the order parameter 
can presumably be calculated by another method. We hope to return to this problem elsewhere. 



5 Special cases 
5.1 BBS 2 model 

Parameters of the general free-fermion model, which correspond to BBS2 model (via the formulas 
(|1.2|) - ()1.9Jl ). are not independent. In particular, in addition to free-fermion condition (jl.lUj) . they 
also satisfy the relation 013024 = 014023- Therefore, one could expect some simplifications of the 
above formulas for tranfer matrix eigenvectors to occur in this case. Furthermore, it is known that 
the eigenvalues of the BBS2 transfer matrix should have polynomial dependence on spectral variable 
t, and that the eigenvectors should not depend on it. In order to verify these properties, let us rewrite 
our formulas in the BBS notation. 

The variables Xp an d Gijip) = 1,2), which were used in the grassmann integral representation 
of the transfer matrix, are expressed in terms of t, x^x , y , y j fjjj fjj cis 

4(t + fifi'xx') 2 + 4(yy' — tfifi') 2 + 8(t + fifi'xx')(yy' — tfifi') cosp 



Xp 

XpGuip) 
X P G 2 2(p) 

X P Gn(p) 



(y + fit)(y> + //) 
8it(y + fifi'x')(y' — ji^'x) sinp 

(y + fity( y > + fi>y 

8it(y — fin'x')(y' + fifi'x) sinp 



(y + ^y' + fi') 2 

16(t 2 + /i 2 /i' 2 (t 2 — x 2 x' 2 ) — y 2 y' 2 ) — 32/i//(t 2 + xx'yy') cosp — 32itfj.fi' (xy + x'y') sinp 

(y + fit) 2 (y> + fu') 2 ' 

In order to write down the eigenvalues of T NS and T R , it is sufficient to express in terms of BBS 
parameters the quantities A^L»i and v v j p p . They are given by (see also 0) 

= ( 1 + ^V L ) U NS (t + t) 
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. d (1 - U L U' L ) i-rfl, 

= y Ly,L n 



v p/p p 



V 

t - 1 



V 



t + t p 



where t p is defined as 



with 



tp 



a p 



a p = 1 — 2/i/i' cosp + fJ?fi' 2 , 
bp = fJ>fjf(xy + x'y') sin p. 

2 /2 , o / / / ,2/22/2 

c p = y y + ^a*/-* 2/2/ cosp + /i /i xx . 

Finally, 'right' and 'left' eigenvectors of T^ 5 and T fl are fully characterized by the functions A ± (p) 
and ^(p), which in the case of the BBS2 model can be written in the following form: 

\±f„\ ^p"F V a p c p ~ bp 



A ± ( P ) 



(y — fifi'x')(y' + fifi'x) isinp' 

A ± ( P ) = 



• d p -\ ^JcipCp 5| 



(y + fifj,'x')(y' — nn'x) i sinp ' 
where 

dp = jj,fj,'(xx' — yy') + (yy 1 — /z 2 // 2 xx') cosp. 
One should note that the functions A + (p), A + (p) do not depend on t, as expected. 

5.2 Ising model 

In the Ising case, another parametrization is typically used. For simplicity, let us consider the 
isotropic model, characterized by the plaquette weight 

W(<7i, a 2 , 0-3, 04) = expji K(a 1 a 2 + o 2 o ? , + (T 3 a 4 + <7 4 <7i)}. 

Parameters of the general free-fermion model, corresponding to this Boltzmann weight, are given by 

cosh 2 K + 1 sinh 2 K 

do = , CI4 = di3 = d24 



cosh 2 K + l' 
sinh K cosh K 



Ol2 — ^23 — ^34 — Ol4 



cosh 2 K + l 

Let us also introduce the function j q , given by the positive root of the equation 

cosh 7g = sinh 2K + sinh" 1 2K — cos q, 
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One can now rewrite the variables Xp an d Gij{p) {hj = 1)2) from the grassmann integral represen- 
tation of T NS and T R in the following way: 



Xp 



sinh 2 2K [(1 + tanh cos p) 2 + sin 2 p] 



XpGxxip) = x P G 2 2(p) 



(cosh 2 K + l) 2 
2i ship sinh 2K (cosh 2K — tanh K cosp) 



X P Gi 2 (p) 



(cosh 2 K + l) 2 
2 sinh 2K 



(cosh 2 X + l) 2 ' 

The eigenvalues and eigenvectors of T NS and T R may be found from 

A-f) = (2 sinh 2i^ 2 exp j \ £* 5(J %j , V* 



e 



"7p 



, ., sinh K cosh T(l + tanh K cosp) 2 + sin 2 pi — e ±7p 

y4 f T) ) — j4 f 77 1 — - - 

' i sin p (cosh 2i^ — tanh K cos p) 

It should be emphasized once again that the transfer matrix eigenvectors, constructed above, auto- 
matically diagonalize the translation operator R as well. Therefore, we believe that these eigenvectors 
may turn out to be useful for the construction of a rigorous proof of the recently obtained formula 
[TU] for Ising spin form factors. 



6 Summary 

We have obtained the transfer matrix eigenvectors of the BBS2 model on a finite lattice, using the 
method of grassmann integration. Our results are exact and explicit, i. e. the eigenvectors are 
expressed in terms of initial lattice variables. Grassmann integral representation for the eigenvectors 
immediately gives their norms and allows to considerably advance in the computation of form factors 
of the BBS2 model. 

The only two things that were actually necessary for our computation are the translational in- 
variance of the model and the representation for the plaquette Boltzmann weight. In this 

respect, the method, developed in the present paper, is quite general and it could be extended to 
free-fermion models with a more complicated configuration space of order parameter, once these are 
found. 
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